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The role of neutron transfers is investigated in the fusion pro- 
cess below the Coulomb barrier by analyzing 32 S+ 90 Zr and 32 S+ 96 Zr 
as benchmark reactions. A full coupled-channel calculation of the 
fusion excitation functions has been performed for both systems by 
using multi-neutron transfer coupling for the more neutron-rich re- 
action. The enhancement of fusion cross sections for 32 S+ 96 Zr is 
well reproduced at sub-barrier energies by NTFus code calculations 
including the coupling of the neutron-transfer channels following 
the Zagrebaev semiclassical model. We found similar effects for 
40 Ca+ 90 Zr and 40 Ca+ 96 Zr fusion excitation functions. 



1. Introduction 

Heavy-ion fusion reactions with colliding neutron-rich nuclei at bom- 
barding energies at the vicinity and below the Coulomb barrier have been 
widely studied [HEIEIIIIEIEIEIEIEI US]- The specific role of multi-step 
neutron-transfers in sub-barrier fusion enhancement still needs to be inves- 
tigated in detail both experimentally EJ [§] and theoretically [HE]. In a 
complete description of the fusion dynamics the transfer channels in stan- 
dard coupled-channel (CC) calculations [H HI El [TUl HI] have to be taken 
into account accurately. It is known that neutron transfers may induce a 
neck region of nuclear matter in-between the interacting nuclei favoring the 
fusion process to occur. 

In low-energy fusion reactions, the very simple one-dimensional barrier- 
penetration model (1D-BPM) is based upon a real potential barrier resulting 
from the attractive nuclear and repulsive Coulomb interactions. For light- 
and medium-mass nuclei, one only assumes that the di-nuclear system fuses 
as soon as it has reached the region inside the barrier i.e. within the potential 
pocket. If the system can evolve with a bombarding energy high enough to 
pass through the barrier and to reach this pocket with a reasonable amount 
of energy, the fusion process will occur after a complete amalgation of the 
colliding nuclei forming the compound nucleus. On the other hand, for 
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sub-barrier energies the di-nuclear system has not enough energy to pass 
through the barrier. In this case, neutron pick-up processes can occur when 
the nuclei are close enough to interact each other significantly E] , if the 
Q-values of neutron transfers are positive. 

It was shown that sequential neutron transfers can lead to the broad dis- 
tributions characteristic of many experimental fusion cross sections. Finite 
Q-value effects can lead to neutron flow and a build up of a neck between 
the target and projectile [3]. The situation of this neck formation of neutron 
matter between the two colliding nuclei could be considered as a "doorway 
state" to fusion. In a basic view, this intermediate state induced a barrier 
lowering. As a consequence, it will favor the fusion process at sub-barrier 
energies and enhance significantly the fusion cross sections. Experimental 
results have already shown such enhancement of the sub-barrier fusion cross 
sections due to the neutron-transfer channels with positive Q-values [2J [5]. 

2. Experimental results 

In order to investigate the role of neutron transfers we further study 32 S 
+ 90 Zr and 32 S + 96 Zr as benchmark reactions. Fig. 1 displays the measured 
fusion cross sections for 32 S + 90 Zr (open circles) and 32 S + 96 Zr (points). 
We present the analysis of excitation functions of evaporation residues (ER) 
cross sections recently measured with high precision (i.e. with small energy 
steps and good statistical accuracy for these reactions [T2]). 

The differential cross sections of quasi-elastic scattering (QEL) at back- 
ward angles were previously measured by the CIAE group [9]. The analy- 
sis of the corresponding BD-QEL barrier distributions (see solid points in 
Fig. 2) already indicated the significant role played by neutron tranfers in 
the fusion processes. 

In Fig. 2 we introduce the experimental fusion-barrier (BD-Fusion) dis- 
tributions (see open poins) obtained for the two reactions by using the 
three-point difference method of Ref. [4j as applied to the data points of 
Ref. [12J plotted in Fig. 1. It is interesting to note that in both cases the 
BD-Fusion and BD-QEL barrier distributions are almost identical up to 
E c . m . « 85 MeV. 

3. Coupled channel analysis 

We have developed a new CC computer code named NTFus [13] by 
taking the neutron transfer channels into account in the framework of the 
semiclassical model of Zagrebaev [6]. We will show that the effect of the 
neutron transfer channels yields a fairly good agreement with the present 
data of sub-barrier fusion cross sections measured for 32 S + 96 Zr, the more 
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Fig. 1. Comparison between the fusion-evaporation (ER) excitation functions of 
32 S+ 90 Zr (open circles) and 32 S+ 96 Zr (points) as a function of the center-of-mass 
energy. The error barrs of the experimental data taken from Ref. |12j represent 
purely statistics uncertainties. 

neutron-rich reaction [12] . This was initially expected from the positive Q- 
values of the neutron transfers as well as from the failure of standard CC 
calculation of quasi-elastic barrier distributions without neutron-transfers 
coupling [9] as shown by the solid line in Fig. 2(b). 

By fitting the present experimental fusion excitation function given in 
Fig. 1 with NTFus CC calculation [13], we will be able to conclude that the 
effect of the neutron transfer channels produces the rather significant en- 
hancement of the sub-barrier fusion cross sections of 32 S + 96 Zr as compared 
to 32 S + 90 Zr. 

A detailed inspection of the 32 S + 90 Zr fusion data presented in Fig. 1 
along with the negative Q-values of their corresponding neutron transfer 
channels lead us to speculate with the absence of a neutron transfer effect 
on the sub-barrier fusion for this reaction. We proceed step by step by 
performing calculations for this reaction with the NTFus code []2] (see 
Fig. 3). 
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Fig. 2. Barrier distributions (BD) from the fusion ER (open circles) cross sections 
[12], plotted in Fig.l, and quasielastic scattering (solid circles) cross sections |9| 
for 32 S+ 90 Zr (a) and 32 S+ 96 Zr (b). The dashed and solid black lines represent un- 
coupled calculations (1D-BPM) and the CC calculations without neutron transfer 
coupling. The red dash-dotted line represents the CC calculations with neutron 
transfer coupling for the 32 S+ 96 Zr reaction. 



4. NTFus CC calculations for 32 S+ 90 ' 96 Zr 

With the semiclassical model developed by Zagrebaev [6] we propose in 
the following discussion to definitively demonstrate the significant role of 
neutron transfers for the 32 S + 96 Zr fusion reaction by fitting its experi- 
mental excitation function with NTFus code [T3] calculations, as shown in 
Fig. 4. The main characteristics of the code are briefly described thereafter. 
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The new oriented object NTFus code [13j . using the Zagrebaev model 
[6] was implemented (at the CIAE) in C++, using the compiler of ROOT 
P3] , following the basic equations of Ref. [15] . 

Let us first remind the values chosen for the deformation parameters and 
the excitation energies that are given in Table 1 [TJ [16j [T7] . The quadrupole 
vibrations of both the 90 Zr and 96 Zr are weak in energy; they lie at compa- 
rable energies. 

Table 1. Excitation energies E^, spins and parities A^ and deformation parameters 

/3 A horn DO US1 [IZI - 
Nucleus E a (MeV) A 71 " (3 X 



32g 


2.230 


2+ 


0.32 




5.006 


3 _ 


0.40 


90 Zr 


2.186 


2+ 


0.09 




2.748 


3" 


0.22 


96 Zr 


1.751 


2+ 


0.08 




1.897 


3 _ 


0.27 



The 96 Zr nucleus presents a complicated situation [18]: its low-energy 
spectrum is dominated by a 2 + state at 1.748 MeV and by a very collec- 
tive [B(E3;3" -> 0+) = 51 W.u.] 3~ state at 1.897 MeV. CC calculations 
explained the larger sub-barrier enhancement as due mainly to the strong 
octupole vibration of the 3~ state in 36 S + 96 Zr [19j . However, the agree- 
ment is not so satisfactory below the barrier for 32 S + 96 Zr (see solid line 
of Fig. 4), as well as for °Ca + 96 Zr [S] and, therefore, there is the need to 
take neutron transfers into account. 

The main functions of the code NTFus are designed to calculate the 
fusion excitation functions with normalized barrier distribution (based on 
experimental data) given by CCFULL [llj . we take the dynamical defor- 
mations into account. 

To take into account the neutron transfers, the NTFus code [13] applies 
the Zagrebaev model [6j to calculate the fusion cross sections af us (E) as a 
function of center-of-mass energy E. Then the fusion excitation function can 
be derived using the following formula [6]: 

Ti(E) = 

/ f( B )TT £ / °<k(E, I, Q) x P HW (B, E + Q, l)dQdB. (1) 

J Wtr ~ J-E 

and 
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(2) 



1=0 



where Ti(E) are the transmission coefficients, E is the energy given 
in the center-of-mass system, B and f(B) are the barrier height and the 
normalized barrier distribution function, ~Phw is the usual Hill- Wheeler 
formula. I is the angular momentum whereas l cr is the critical angular 
momentum as calculated by assuming no coupling (well above the barrier). 
ctk(E,l,Q) and Qo(k) are, respectively, the probabilities and the Q-values 
for the transfers of k neutrons. And 1/Nt r is the normalization of the total 
probability taking into account the neutron transfers. 

The NTFus code [13] uses the ion-ion potential between two deformed 
nuclei as developped by Zagrebaev and Samarin in Ref. [15]. Either the 
standard Woods-Saxon form of the nuclear potential or a proximity poten- 
tial [20] can be chosen. The code is also able to predict fusion cross sections 
for reactions induced by halo projectiles; for instance 6 He + 64 Zn [21]. In 
the following, only comparisons for 32 S + 90 Zr and 32 S + 96 Zr are discussed. 

For the high-energy part of the 32 S + 90 Zr excitation function, one can 
notice a small over-estimation of the fusion cross sections at energies above 
the barrier up to the point used to calculate the critical angular momentum. 
This behavior can be observed at rather high incident energies - i.e. between 
about 82 MeV and 90 MeV (shown in Fig. 3 for 32 S + 90 Zr reaction). We 
want to stress that the corrections do not affect our conclusions that the 
transfer channels have a predominant role below the barrier for 32 S + 96 Zr 
reaction, as shown in Fig. 4. 

As expected, we obtain a good agreement with calculations not tak- 
ing any neutron transfer coupling into account for 32 S + 90 Zr as shown by 
the solid line of Fig. 3 (the dashed line are the results of calculations per- 
formed without any coupling). On the other hand, there is no significant 
over-estimation at sub-barrier energies. As a consequence, it is possible to 
observe the strong effect of neutron transfers on the fusion for the 32 S + 96 Zr 
reaction at sub-barrier energies. Moreover, the barrier distribution function 
f(B) extracted from the data contains the information of the neutron trans- 
fers. These information are also contained in the transmission coefficients, 
which are the most important parameters for the fusion cross sections to be 
calculated accurately. The f(B) function as calculated with the three-point 
formula [3j will mimic the differences induced by the neutron transfer tak- 
ing place in sub-barrier energies where the cross section variations are very 
small (only visible if a logarithm scale is employed for the fusion excitation 
function). It is interesting to note that the Zagrebaev model [6] implies a 
modification of the Hill- Wheeler probability and does not concern the bar- 
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Fig. 3. Fusion-evaporation (ER) excitation function for 32 S+ 90 Zr. The solid points 
are the experimental data [12] (see Fig.l). The dashed and solid lines are the un- 
coupled calculations, and CC calculations without neutron transfers, respectively. 
The arrow indicates the position of the Coulomb barrier for 32 S+ 90 Zr as given by 
the 1D-BPM model (see Fig. 2). 




rier distribution function f(B). Finally, the code allows us to perform each 
calculation by taking the neutron transfers into account or not. 

Table 2. Q-values in MeV for neutron pickup transfer channels from ground state 
to ground state for 32 S+ 90 Zr and 32 S+ 96 Zr, respectevly. 
System +ln +2n +3n +4n 



3i S+ 90 Zr -3.33 -1.229 -6.59 -6.319 
32 S+ 96 Zr 0.788 5.737 4.508 7.655 

The calculation with the neutron transfer effect is performed here up to 
the channel +4n (k=4), but we have seen that we obtain the same overall 
agreement with data up to channels +5n and +6n. The Q-values and the 
separation energies for the 96 Zr nucleus used for this calculation (solid line 
in Fig. 4) are displayed in the Tables 2 and 3, respectively. 
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Fig. 4. Fusion-evaporation (ER) excitation function for 32 S+ 96 Zr. The solid points 
are the experimental data |T2] (see Fig.l). The dashed, solid, and dotted lines 
are the uncoupled calculations, and CC calculations without and with neutron 
transfers, respectively. The arrow indicates the position of the Coulomb barrier for 
32 S+ 96 Zr as given by the id-BPM model (see Fig. 2). 

Table 3. Separation energies in MeV of each neutron for 96 Zr. 
I s * neutron 2 nd neutron 3 rd neutron 4 th neutron 



7.854 6.463 8.230 6.733 



As we can see on Fig. 4, the solid line representing standard CC cal- 
culations without the neutron transfer coupling (the dotted line is given 
for uncoupled calculations) does not fit the experimental data well at sub- 
barrier energies. On the other hand, the dotted line displaying NTFus 
calculations taking the neutron transfer coupling into account agrees per- 
fectly well with the data. As expected, the Zagrebaev semiclassical model's 
correction applied at sub-barrier energies enhances the calcutated cross sec- 
tions. Moreover, it allows to fit the data reasonably well and therefore 
illustrates the strong effect of neutron transfers for the fusion of 32 S + 96 Zr 
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at subbarrier energies. 

The present full CC analysis of 32 S + 96 Zr fusion data [12j using NTFus 
|13j confirms perfectly well first previous CC calculations [B] describing well 
the earlier 40 Ca + 90 ' 96 Zr fusion data [S] and, secondly, very recent fragment- 
7 coincidences measured for 40 Ca + 96 Zr multi-neutron transfer channels 

US]. 

5. Summary and conclusions 

We have investigated the fusion process (excitation functions and ex- 
tracted barrier distributions [12]) at near- and sub-barrier energies for the 
two neighbouring reactions 32 S + 90 Zr and 32 S + 96 Zr. For this purpose a 
new computer code named NTFus [13] has been developped by taking the 
coupling of the multi-neutron transfer channels into account by using the 
semiclassical model of Zagrebaev [B]. 

The effect of neutron couplings provides a fair agreement with the present 
data of sub-barrier fusion for 32 S + 96 Zr. This was initially expected from 
the positive Q-values of the neutron transfers as well as from the failure of 
previous CC calculation of quasi-elastic barrier distributions without cou- 
pling of the neutron transfers [9J. With the agreement obtained by fitting 
the present experimental fusion excitation function and the CC calculation 
at sub-barrier energies, we conclude that the effect of the neutron transfers 
produces a rather significant enhancement of the sub-barrier fusion cross 
sections of 32 S + 96 Zr as compared to 32 S + 90 Zr. At this point we did 
not try to reproduce the details of the fine structures observed in the fu- 
sion barrier distributions. We believe that to achieve this final goal it will 
first be necessay to measure the neutron transfer cross sections to provide 
more information on the coupling strenght of neutron transfer because its 
connection with fusion is not yet fully understood [18) . 
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